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Abstract 

We have analytically obtained all the density matrix elements up to six lattice 
sites for the spin-1/2 Heisenberg XXZ chain at A = 1/2. We use the multiple inte- 
gral formula of the correlation function for the massless XXZ chain derived by Jimbo 
and Miwa. As for the spin-spin correlation functions, we have newly obtained the 
fourth- and fifth- neighbour transverse correlation functions. We have calculated all 
the eigenvalues of the density matrix and analyze the eigenvalue-distribution. Using 
these results the exact values of the entanglement entropy for the reduced density ma- 
trix up six lattice sites have been obtained. We observe that our exact results agree 
quite well with the asymptotic formula predicted by the conformal field theory. 
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1 Introduction 



The spin- 1/2 antiferromagnetic Heisenberg XXZ chain is one of the most fundamental models 
for one-dimensional quantum magnetism, which is given by the Hamiltonian 

j=—ao 

where S" = o~J/2 with a? being the Pauli matrices acting on the j-th site and A is the 
anisotropy parameter. For A > 1, it is called the massive XXZ model where the system is 
gapful. Meanwhile for — 1 < A < 1 case, the system is gap less and called the massless XXZ 
model. Especially we call it XXX model for the isotropic case A = 1. 

The exact eigenvalues and eigenvectors of this model can be obtained by the Bethe Ansatz 
method [TJ [2]. Many physical quantities in the thermodynamic limit such as specific heat, 
magnetic susceptibility, elementary excitations, etc., can be exactly evaluated even at finite 
temperature by the Bethe ansatz method [2]. 

The exact calculation of the correlation functions, however, is still a difficult problem. The 
exceptional case is A = 0, where the system reduces to a lattice free-fermion model by the 
Jordan- Wigner transformation. In this case, we can calculate arbitrary correlation functions 
by means of Wick's theorem [3111]. Recently, however, there have been rapid developments in 
the exact evaluations of correlation functions for A / case also, since Kyoto Group (Jimbo, 
Miki, Miwa, Nakayashiki) derived a multiple integral representation for arbitrary correlation 
functions. Using the representation theory of the quantum affine algebra U q {sl2), they first 
derived a multiple integral representation for massive XXZ antiferromagnetic chain in 1992 
[5J E], which is before long extended to the XXX case [H [5] and the massless XXZ case 
[9]. Later the same integral representations were reproduced by Kitanine, Maillet, Terras 
[TO] in the framework of Quantum Inverse Scattering Method. They have also succeeded in 
generalizing the integral representations to the XXZ model with an external magnetic field 
[10] . More recently the multiple integral formulas were extended to dynamical correlation 
functions as well as finite temperature correlation functions [TTj [T2"l [T5| [T3] . In this way 
it has been established now the correlation functions for XXZ model are represented by 
multiple integrals in general. However, these multiple integrals are difficult to evaluate both 
numerically and analytically. 

For general anisotropy A, it has been shown that the multiple inetegrals up to four- 
dimension can be reduced to one- dimensional integrals [T51 [TBI El HSl [333 EDI [21]. As a 
result all the density matrix elements within four lattice sites have been obtained for general 
anisotropy [21] . To reduce the multiple integrals into one-dimension, however, involves hard 
calculation, which makes difficult to obtain correlation functions on more than four lattice 
sites. On the other hand, at the isotropic point A = 1, an algebraic method based on 
qKZ equation has been devised [22] and all the density matrix elements up to six lattice 
sites have been obtained [2"5| |2"3]. Moreover, as for the spin-spin correlation functions, up to 
seventh-neighbour correlation (SfS§) for XXX chain have been obtained from the generating 
functional approach [251 [26] . It is desirable that this algebraic method will be generalized to 
the case with A^l. Actually, Boos, Jimbo, Miwa, Smirnov and Takeyama have derived an 
exponential formula for the density matrix elements of XXZ model, which does not contain 
multiple integrals [23 [2H1 [2S1 EDI EE]- It, however, seems still hard to evaluate the formula 
for general density matrix elements. 
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Among the general A^O, there is a special point A = 1/2, where some intriguing prop- 
erties have been observed. Let us define a correlation function called Emptiness Formation 
Probability (EFP) [8j which signifies the probability to find a ferromagnetic string of length 
n: 




The explicit general formula for P(n) at A = 1/2 was conjectured in [33] 

k=0 v ' 

which is proportional to the number of alternating sign matrix of size n x n. Later this 
conjecture was proved by the explicit evaluation of the multiple integral representing the 
EFP [M]. Remarkably, one can also obtain the exact asymptotic behavior as n — > oo 
from this formula, which is the unique valuable example except for the free fermion point 
A = 0. Note also that as for the longitudinal two-point correlation functions at A = 1/2, 
up to eighth-neighbour correlation function (SfSg) have been obtained in [32] by use of the 
multiple integral representation for the generating function. Most outstanding is that all the 
results are represented by single rational numbers. These results motivated us to calculate 
other correlation functions at A = 1/2. Actually we have obtained all the density matrix 
elements up to six lattice sites by the direct evaluation of the multiple integrals. All the 
results can be written by single rational numbers as expected. A direct evaluation of the 
multiple integrals is possible due to the particularity of the case for A = 1/2 as is explained 
below. 



2 Analytical evaluation of multiple integral 

Here we shall describe how we analytically obtain the density matrix elements at A = 1/2 
from the multiple integral formula. Any correlation function can be expressed as a sum of 
density matrix elements Pel'--- 'en, which are defined by the ground state expectation value of 
the product of elementary matrices: 

pt::i: = (Ef ei ---E^) 1 (2.1) 



where E? are 2x2 elementary matrices acting on the j-th site as 



Et+ ~ (o o) M " \ + s *> *T " (o i) w -\~^ 
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The multiple integral formula of the density matrix element for the massless XXZ chain 
reads [H] 

pfv->< = ( v \-n(n-i)/2 f°° dx! /*°° dx^ T t sinh(x a - x b ) 

J-oo 2vr 2ir l± b smh[(x a -x b -if ab ir)i>] 

x " sinh^- 1 [(xfc + ztt/2>] sinh"-^ [fo - m/2)u] 
f-*- cosh™Xfc 

k=l 

where the parameter v is related to the anisotropy as A = cos7rz/ and f ab and y^ are 
determined as 

f ab = (1 + sign[(s' - a + l/2)(s' - b + l/2)])/2, 

2/1 > vi > ■ ■ ■ > y s >, 4 = + 

y s ' + i > ■■■ > y n , e n+ i_ yj = -. (2.3) 

In the case of A = 1/2, namely v = 1/3, the significant simplification occurs in the multiple 
integrals due to the trigonometric identity 

sinh(x a — x b ) = 4sinh[(x a — x b )/3) sinh[(a; a — x b +in)/3] smh[(x a —x b —i7r)/3]. (2.4) 

Actually if we note that the parameter f ab takes the value or 1, the first factor in the 
multiple integral at v — 1/3 can be decomposed as 

sinh(x a - x b ) _ ^ ^ h / x a - x b \ ( x a - x b + in 



sinh[(x a — x b — in)/ 3] \ 3 / \ 3 

= -1 + ue 2 3 (Xa - Xb) + u- x e-^ Xa - Xb \ (2.5) 

sinh(x a - x b ) ( x a - x b + in\ . f x a - x b - in 
4 sinn smh 



sinh[(x a — x b )/3] \ 3 

= I + e U x °- x i>) + e ~*( Xa ~ Xb \ (2.6) 

where to = e t7T ^ 3 . Expanding the trigonometoric functions in the second factor into exponen- 
tials 

sinh^ 1 [(x + in/2)/3] sinh n " y [(x - in/ 2)/ 3] 

= 2 l - n (oo^e^ - W -V2 e -*/3jW-l ( w -l/2 e */3 _ w l/2 e - 3 ,/3jn-» 

Z I \ m 



2±-n ^2( — l) l+m ^ ( n ^) [J 2/-i+m-(n+l)/2 e |(n-2Z-2m-l)a: 

i=0 m=0 ^ 



we can explicitly evaluate the multiple integral by use of the formula 

00 e ax dx nn i „ / n + a n — a\ _ , . . „. 

2"~ 1 J B , , Re(n±ct)>0, (2.8) 



oo 



cosh™ x 

where B(p, q) is the beta function defined by 



B(p,q)= [ t p -\l-t) q - l dt, Re(p),Re(g) > 0. 



(2.9) 
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Table 1: Comparison with the asymptotic formula of the transverse correlation function 









(SfSl) 


(SfSI) 


(S?S$) 


Exact 


-0.156250 


0.0800781 


-0.0671234 


0.0521997 


-0.0467664 


Asymptotics 


-0.159522 


0.0787307 


-0.0667821 


0.0519121 


-0.0466083 



In this way we have succeeded in calculating all the density matrix elements up to six lattice 
sites. All the results are represented by single rational numbers, which are presented in 
Appendix A. As for the spin-spin correlation functions, we have newly obtained the fourth- 
and fifth-neighbour transverse two-point correlation function 

5 



(5?5f> 



-0.15625, 



0.080078125, 



32 
41 

512 " 
4399 

~ 65536 " 
1751531 

33554432 
3213760345 



-0.0671234130859375, 
= 0.0521996915340423583984375, 

-0.046766368104727007448673248291015625. 



68719476736 

The asymptotic formula of the transverse two-point correlation function for the massless 
XXZ chain is established in 1351 [361 



(S?S? + „> ~ A X ( V ) 



ell 

VP 



-Mrf) 



1 



5(1 -nf 



i 

2^(1 - V) 



2-2r] 



V+ 
Tl V 

-I V 



T + 



7} 



I -is, 



20Fr 
r 



2— 2r? / _ 



exp 



smh(nt) 



2-2r) 



-i *H 



2v^Fr 



2-2n J _ 



exp 



sinh(t) cosh[(l — r))t] 

cosh(2r]t)e~ 2t 



-2t 



dt 

T 



2 sinh(?7t) sinh(t) cosh[(l — rj)t\ 



-2t 



dt 
T 



. ... (2.10) 

smh[T]t) rj J t 

which produces a good numerical value even for small n as is shown in Table [U Note that the 
longitudinal correlation function was obtained up to eighth-neighbour correlaion (SfSg) from 
the multiple integral representation for the generating function [32J. Note also that up to 
third-neighbour both longitudinal and transverse correlation functions for general anisotropy 
A were obtained in [21j . 



3 Reduced density matrix and entanglement entropy 

Below let us discuss the reduced density matrix for a sub-chain and the entanglement entropy. 
The density matrix for the infinite system at zero temperature has the form 

p T = |GS)(GS|, (3.1) 
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Figure 1: Eigenvalue-distribution of density matrices 



Table 2: Entanglement entropy S(n) of a finite sub-chain of length n 



5(1) 


5(2) 


5(3) 


5(4) 


1 


1.3716407621868583 


1.5766810784924767 


1.7179079372711414 



5(5) 



5(6) 



1.8262818282012363 1.9144714710902746 



where |GS) denotes the ground state of the total system. We consider a finite sub-chain of 
length n, the rest of which is regarded as an environment. We define the reduced density 
matrix for this sub-chain by tracing out the environment from the infinite chain 



Pn 



tr E PT 



(3.2) 



We have numerically evaluate all the eigenvalues uj a (a — 1, 2, • • • , 2 n ) of the reduced density 
matrix p n up to n = 6. We show the distribution of the eigenvalues in Figure [H The 
distribution is less degenerate comapared with the isotropic case A = 1 shown in [24] . In the 
odd n case, all the eigenvalues are two-fold degenerate due to the spin-reverse symmetry. 

Subsequently we exactly evaluate the von Neumann entropy (Entanglement entropy) 
defined as 



5 (n) = -trp n log 2 p n = -^uj a log 2 u a . 



(3.3) 



a=l 



The exact numerical values of S(n) up to n = 6 are shown in Table [2j By analyzing the 
behaviour of the entanglement 5(n) for large n, we can see how long quantum correlations 
reach [3?]. In the massive region A > 1, the entanglement entropy will be saturated as n 
grows due to the finite correlation length. This means the ground state is well approximated 
by a subsystem of a finite length corresponding to the large eigenvalues of reduced density 
matrix. On the other hand, in the massless case — 1 < A < 1, the conformal field theory 
predict that the entanglement entropy shows a logarithmic divergence [38] 



5(n) ~ -log 2 n + k A . 



(3.4) 
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Figure 2: Entanglement entropy Sin) of a finite sub-chain of length n 



Our exact results up to n 



[2j We estimate the numerical value of the constant term &a=i/2 as &a=i/2 ~ S(6) — | log 2 6 



6 agree quite well with the asymptotic formula as shown in Figure 

as ki 

1.0528. This numerical value is slightly smaller than the isotropic case A = 1, where the 
constant &a=i is estimated as /ca=i ~ 1.0607 from the exact data for S(n) up to n = 6 [24J. 
At free fermion point A = 0, the exact asymptotic formula has been obtained in [39] 



S(n) ~ -log 2 n + k A =o, 



"A=0 



1/3 



dt 



e-* ■ 1 coah(t/2) \ 

3t tsinh 2 (t/2) 2sinh 3 (t/2) J 7 



In this case the numerical value for the constant term is given by &a= 



(3.5) 
1.0474932144- ■ •. 



4 Summary and discussion 

We have succeeded in obtaining all the density matrix elements on six lattice sites for XXZ 
chain at A = 1/2. Especially we have newly obtained the fourth- and fifth-neighbour 
transverse spin-spin correlation functions. Our exact results for the transverse correlations 
show good agreement with the asymptotic formula established in [SHIES]- Subsequently we 
have calculated all the eigenvalues of the reduced density matrix p n up to n = 6. From these 
results we have exactly evaluated the entanglement entropy, which shows a good agreement 
with the asymptotic formula derived via the conformal field theory. Finally, we remark 
that similar procedures to evaluate the multiple integrals are also possible at v = 1/n for 
n = 4, 5, 6, • • • , since there are similar trigonometric identities as (12. 4p . We will report the 
calculation of correlation functions for these cases in subsequent papers. 
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Appendix A Density matrix elements up to n = 6 

In this appendix we present all the independent density matrix elements defined in eq. (12.11) 
up to n = 6. Other elements can be computed from the relations 



p e l>- > e 



if E f ^E^ 



e' ■•• e' 



+ p: 



— ei)— ,-£n. 

I p e l'"' ' € ro 

, _ 61,— ,fn; 



and the formula for the EFP [33l El] 

/'i//) Pic:' 



n— 1 



2-n 



fc=0 



(3fc + l)! 
(n + fe)! ' 



(A.l) 

(A.2) 
(A.3) 



(A.4) 



Appendix A.l n < 4 



^+4-4 

r +4~ 
p— H 



16 
221 

"8192 " 
289 

"32768 
2005 

" 32768 
1393 

65536 " 
3091 



-0.3125, 
= -0.0269775, 
= -0.00881958, 
= -0.0611877, 
0.0212555, 



r 32768 
Appendix A.2 n = 5 



0.0943298. 



p-++ 
p-+++ 

r +4— 4- 

p4— H- 
r + -4-4- 

p— ++ 

r + -+- 

p-4~+ 

-t _l I 



41 



0.0800781, 



512 
1579 

65536 =°- 0240936 ' 
1037 

= 0.0632935, 

16384 



3821 

" 65536 
4883 

32768 ~ 



-0.0583038, 



0.149017, 



pi 



-+++ 



14721 

"8388608 " 
48987 

"33554432 
179699 

33554432 " 
165155 

33554432 
31069 



-0.00175488, 



-0.00145993, 



0.00535545, 
0.004922, 



0.0148149, 



2097152 

196569 

= 0.0117164, 

16777216 



p-+++- 
r ++-+-\ 

p-+++- 
r ++++- 

P+-++- 

J I I ! 



P++-+H 
4-4— 1-4 



p +++ 

r +-+-+ 

p h++ 

*++-+- 



37335 

16777216 " 
13911 

33554432 " 
120337 

"16777216 
168313 

16777216 " 
411583 

"16777216 
281271 

"33554432 



0.00222534, 
0.00041458, 



-0.00717264, 



0.0100322, 



-0.0245323, 
-0.00838253, 
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p, 

p: 



-+-++ 
-++- 



p 



+-++- 



p-+++- 
r +-+-+ 



79673 

33554432 " 
1261655 

"33554432 
1575515 

33554432 " 
1366619 



0.00237444, 
= -0.0376002, 



0.046954, 



33554432 
Appendix A. 3 n = 6 



0.0407284. 



p-+-++ 
r ++-+- 



-++ 



1441787 

"33554432 
59459 

2097152 " 
696151 

"33554432 



= -0.0429686, 
0.0283523, 
= -0.0207469, 



P-+++++ 
r +-++++ 



P-+++++ 
r +++-++ 



1546981 



++-+++ 



P++-+++ 
r +++-++ 



-++ 



p — ++++ 
- r +-+++- 



p — ++++ 
r ++-+-+ 



p — ++++ 
r +++ — + 



-++++ 



P-+-+++ 
-*+-+-++ 



P-+-+++ 



P-+-+++ 
^+++ — + 



-+++ 



P-++-++ 



P-++-++ 
r ++ — ++ 



P-++-++ 



34359738368 
2366275 

"34359738368 
284577 

"34359738368 
20086627 

"68719476736 
10295153 

"68719476736 
35087523 

"68719476736 
214080091 

"68719476736 
57522267 

"68719476736 
154538459 

"68719476736 
6708473 

8589934592 " 
3860673 

34359738368 
12211375 



-0.0000450231, 
-0.0000688677, 
-8.28228 x 10~ 6 , 
-0.000292299, 
-0.000149814, 
-0.000510591, 
-0.00311528, 
-0.000837059, 
= -0.00224883, 
0.000780969, 



P-+++++ 



P-+++++ 



1073741824 
56183761 

17179869184 
606065059 

68719476736 " 
108202041 

"34359738368 
7860495 

1073741824 " 
1044016671 

68719476736 " 
676957849 

68719476736 " 
6581795 

1073741824 



= 0.00011236, 
0.0113727, 



= 0.00327033, 
= 0.00881941, 
= -0.00314909, 
0.00732066, 
= 0.0151924, 
= 0.00985103, 
0.00612977, 



Pi 



P+-++++ 
r +++-++ 



P 



++- 



+ — +++ 



+-++- 



p — ++++ 
r ++ — ++ 



p — ++++ 



p — ++++ 
r +++-+- 



p: 



+-+++ 
— +++ 

f-+++ 
-++-+ 



p-+-+++ 
r ++ — ++ 



p-+-+++ 
1 ++-++- 



p-+-+++ 
^+++-+- 



p-++-++ 
r +-+-++ 



p-++-++ 
r +-+++- 



p-++-++ 
r ++-+-+ 



p-++-++ 
r +++ — + 



5095899 

68719476736 " 
2455833 

68719476736 " 
2927709 

17179869184 " 
19268565 

68719476736 " 
17781349 

34359738368 " 
48421023 

34359738368 " 
88171589 

34359738368 " 
56776545 

34359738368 " 
60809571 

68719476736 " 
33366621 

"68719476736 
85706851 

~ 17179869184 
332557469 

"34359738368 
430452959 

"68719476736 
123612511 

"34359738368 
70061315 

34359738368 " 
591759525 

"34359738368 
367905053 

"68719476736 
988973861 

"68719476736 
363618785 

68719476736 " 



= 0.0000741551, 
= 0.0000357371, 
= 0.000170415, 
= 0.000280395, 
= 0.000517505, 
= 0.00140924, 
= 0.00256613, 
= 0.00165241, 
= 0.000884896, 
= -0.000485548, 
= -0.0049888, 
= -0.0096787, 
= -0.00626391, 
= -0.0035976, 
= 0.00203905, 
= -0.0172225, 
= -0.00535372, 
= -0.0143915, 
= 0.00529135, 
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185522333 

~ 34359738368 
103539423 

"8589934592 = 
267901987 

"34359738368 
309855965 

"68719476736 
296882461 

"68719476736 
92176287 

34359738368 " 
972245985 

"68719476736 
211696415 

"17179869184 
1196499417 

34359738368 " 
1108384987 



34359738368 
347202525 

17179869184 
46285135 

8589934592 " 
19939391 

17179869184 " 
1018463205 

68719476736 " 
277721503 

34359738368 ~ 
369408975 

"17179869184 
880560357 

"68719476736 
113631201 

17179869184 " 
548645951 

34359738368 ~ 
1719255909 

34359738368 " 
1565770597 

34359738368 " 
2117554719 

"68719476736 



= -0.00539941, 
= -0.0120536, 
= -0.00779697, 
= -0.004509, 
= -0.00432021, 
: 0.00268268, 
= -0.014148, 
= -0.0123224, 
0.0348227, 
0.0322582, 



= 0.0202098, 



0.0053883, 
= 0.00116063, 
= 0.0148206, 
= 0.00808276, 
= -0.0215024, 
= -0.0128138, 
0.00661421, 
0.0159677, 
0.0500369, 
0.0455699, 



= -0.0308145. 



PI 



P-+++-+ 
- r +-+++- 



++- 



+ — +++ 



P 



+-++- 



p+ — +++ 
r + — +++ 



p+ — +++ 
r +-++-+ 



p+ — +++ 
r +++ — + 



P+-+-++ 
r +-++-+ 



P\ 

p: 



--++-+ 
-++-+ 

— +++ 
— ++- 

— +++ 
-++ — 



p — +++ 
r +++ — 



-++ 



p — +-++ 
r ++-+ — 



p — ++-+ 
r +-+-+- 



p — ++—\ 
r ++ — +- 



-+++- 



p — +++- 
r ++ — + 



+-+-+- 



901633567 

68719476736 
38524625 

8589934592 = 
12750645 

1073741824 " 
29410257 

17179869184 
35985105 

8589934592 " 
202646807 

34359738368 " 
217687057 

17179869184 " 
78922695 

17179869184 " 
2209522727 

"68719476736 
530683585 

17179869184 ~ 
268623007 

"68719476736 
136974885 

"68719476736 
18442085 

"68719476736 
1454513249 

"68719476736 
335265249 



= 0.0131205, 



68719476736 
1104236607 

34359738368 " 
876924641 

"68719476736 
292857807 

"17179869184 
377925345 

"68719476736 
5350158879 

"68719476736 
3059753503 

"68719476736 



= 0.00448486, 
= 0.011875, 
= 0.0017119, 
= 0.00418922, 
= 0.0058978, 
= 0.0126711, 
= 0.00459391, 
= -0.0321528, 
0.0308899, 
= -0.00390898, 
= -0.00199325, 
= -0.000268368, 
= -0.021166, 
= -0.00487875, 



= 0.0321375, 
= -0.0127609, 
= -0.0170466, 
= -0.00549954, 
= -0.0778551, 
= -0.0445253, 



10 



References 



[i 

[2 

[3] 
[4] 
[5 

[e; 

[7 



H.A. Bethe, Z. Phys. 71 (1931) 205. 

M. Takahashi, Thermodynamics of One- Dimensional Solvable Models, Cambridge Uni- 
versity Press, Cambridge, 1999. 

E. Lieb, T. Schultz, D. Mattis, Ann. Phys. (N.Y.) 16 (1961) 407. 

B.M. McCoy, Phys. Rev. 173 (1968) 531. 

M. Jimbo, K. Miki, T. Miwa, A. Nakayashiki, Phys. Lett. A 168 (1992) 256. 

M. Jimbo, T. Miwa, Algebraic Analysis of Solvable Lattice Models, CBMS Regional Con- 
ference Series in Mathematics vol.85, American Mathematical Society, Providence, 1994. 



A. Nakayashiki, Int. J. Mod. Phys. A 9 (1994) 5673. 

V.E. Korepin, A. Izergin, F.H.L. Essler, D. Uglov, Phys. Lett. A 190 (1994) 182. 

[9] M. Jimbo, T. Miwa, J. Phys. A: Math. Gen. 29 (1996) 2923. 

[10] N. Kitanine, J.M. Maillet, V. Terras, Nucl. Phys. B 567 (2000), 554. 

[11] N. Kitanine, J.M. Maillet, N.A. Slavnov, V. Terras, Nucl. Phys. B 729 (2005) 558. 

[12] F.G6hmann, A. Kliimper, A. Seel, J. Phys. A: Math. Gen 37 (2004) 7625. 

[13] F.G6hmann, A. Kliimper, A. Seel, J. Phys. A: Math. Gen 38 (2005) 1833. 

[14] K. Sakai, "Dynamical correlation functions of the XXZ model at finite temperature", 
|cond-mat/0703319[ 

[15] H.E. Boos, V.E. Korepin, J. Phys. A: Math. Gen. 34 (2001) 5311. 

[16] H.E. Boos, V.E. Korepin, "Evaluation of integrals representing correlators in XXX 
Heisenberg spin chain" in. MathPhys Odyssey 2001, Birkhauser, Basel, (2001) 65. 

[17] H.E. Boos, V.E. Korepin, Y. Nishiyama, M. Shiroishi, J. Phys. A: Math. Gen 35 (2002) 
4443. 

[18] K. Sakai, M. Shiroishi, Y. Nishiyama, M. Takahashi, Phys. Rev. E 67 (2003) 065101. 

[19] G. Kato, M. Shiroishi, M. Takahashi, K. Sakai, J. Phys. A: Math. Gen. 36 (2003) L337. 

[20] M. Takahashi, G. Kato, M. Shiroishi, J. Phys. Soc. Jpn, 73 (2004) 245. 

[21] G. Kato, M. Shiroishi, M. Takahashi, K. Sakai, J. Phys. A: Math. Gen. 37 (2004) 5097. 

[22] H.E. Boos, V.E. Korepin, F.A. Smirnov, Nucl. Phys. B 658 (2003) 417. 

[23] H.E. Boos, M. Shiroishi, M. Takahashi, Nucl. Phys. B 712 (2005) 573. 

[24] J. Sato, M. Shiroishi, M. Takahashi, J. Stat. Mech. 0612 (2006) P017. 



11 



[25] J. Sato, M. Shiroishi, J. Phys. A: Math. Gen. 38 (2005) L405. 



[26] J. Sato, M. Shiroishi, M. Takahashi, Nucl. Phys. B 729 (2005) 441, |hep-th/ 0507290. 

[27] H.E. Boos, M. Jimbo, T. Miwa, F. Smirnov, Y. Takeyama, Algebra Anal. 17 (2005) 
115. 

[28] H.E. Boos, M. Jimbo, T. Miwa, F. Smirnov, Y. Takeyama, Commun. Math. Phys. 261 
(2006) 245. 

[29] H.E. Boos, M. Jimbo, T. Miwa, F. Smirnov, Y. Takeyama, J. Phys. A: Math. Gen. 38 

(2005) 7629. 

[30] H.E. Boos, M. Jimbo, T. Miwa, F. Smirnov, Y. Takeyama, Lett. Math. Phys. 75 (2006) 
201. 

[31] H.E. Boos, M. Jimbo, T. Miwa, F. Smirnov, Y. Takeyama, Annales Henri Poincare 7 

(2006) 1395. 

[32] N. Kitanine, J.M. Maillet, N.A. Slavnov, V. Terras, J. Stat. Mech. 0509 (2005) L002. 

[33] A.V. Razumov, Yu.G. Stroganov, J. Phys. A: Math. Gen. 34 (2001) 3185. 

[34] N. Kitanine, J.M. Maillet, N.A. Slavnov, V. Terras, J. Phys. A: Math. Gen. 35 (2002) 
L385. 

[35] S. Lukyanov, A. Zamolodchikov, Nucl. Phys. B 493 (1997) 571. 
[36] S. Lukyanov, V. Terras, Nucl. Phys. B 654 (2003) 323. 

[37] G. Vidal, J.I. Latorre, E. Rico, A. Kitaev, Phys. Rev. Lett. 90 (2003) 227902. 
[38] C. Holzhey, F. Larsen, F. Wilczek, Nucl. Phys. B 424 (1994) 443. 
[39] B.-Q. Jin, V.E. Korepin, J. Stat. Phys. 116 (2004) 79. 



12 



